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FOREWORD 


This  final  report  presents  the  results  and  conclusions  of  the  Study  of 
Infrared  Emission  in  Heated  Air  carried  out  by  the  Aerodynamic  Research 
Department  of  Cornell  Aeronautical  Laboratory.  The  research  was  supported 
under  Army  Ordnance  Prime  Contract  DA- 1 1  -  022 -O  RD- 3  1 30  by  Subcontract 
Purchase  Order  C  303681  G,  from  the  Bendix  Aviation  Corporation. 
It  IS  the  fourth  report  on  this  program,  having  been  preceded  by  CAL  Report 
Numbers  QM-  1373-A-1.  QM-  1373-A-2,  and  QM-1373-A-3.  with  the  same 
title.  It  IS  understood  that  this  research  is  to  be  continued  under  the  direct 
sponsorship  of  ARGMA. 

The  relationship  of  the  radiation  measurements  contained  herein  to 
re-  entry  has  been  reported  in  the  Proceedings  of  the  recent  AMRAC  meeting 
in  Denver,  Colorado,  May  1961. 
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ABSTRACT 


This  sepost  presents  the  results  of  a  research  progran,  whose  objective 
ts  to  formulate  a  comprehensive  picture  of  the  radiation  from  re-entry  bodies. 
The  program  is  divided  into  two  phases.  the  experimental  determination  of 
the  radiative  properties  of  the  opticaliy-active  gaseous  species,  and  the 
analysis  of  the  blunt-body  flow  field,  in  which  the  f.nite  rate  chemistry  of  the 
air  is  coupled  with  the  aerodynamics  of  the  flow. 

Two  particular  radiating  systems  were  studied,  and  the  results  are  re¬ 
ported  herein.  Measurements  have  been  ccmoleted  on  the  continuum  radiation 

k 


from  the  recombination  reaction  NO  +  O 


NO 2  +  h-v>  in  the  wavelength 


range  0.  5  to  1.  0  microns.  The  value  of  the  rate  constant  at  3750-K,  k  = 

4.4X10  cm  mole  ‘  ‘ ,  is  some  20  times  lower  than  that  measured  at 

room  temperature.  The  spectral  intensity  of  this  radiation  is  presented. 

The  spectrum  of  the  first  positive  band  system  of  nitrogen  between  0.  7 
.nd  1.  5  microns  has  been  obtained.  Absolute  intensities  are  presented  for  a 
range  of  tem.peratu res  and  densities.  The  temperature-dependence  of  the 
radiation  is  shown  to  be  in  good  agreement  with  predicted  values.  The  results 
of  this  work  are  compared  with  those  of  Keck,  et  al  (Annals  of  Physics  7,  1, 
(1959)  )  which  were  based  on  measurements  made  in  air.  The  measured  inten¬ 
sities  from  pure  nitrogen  are  shown  to  be  less  by  a  factor  of  five. 

The  solution  for  the  inviscid.  nonequilibr.um  flow  around  blunt  axisym- 
bodies  (i.  e.  the  bow  shocK  solution)  has  been  programmed  for  an  IBM- 
704  computer.  In  addition,  a  similar  finite  ■  rate  normal-shock  solution  has  been 
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programmed  to  study  the  overall  characteristics  of  the  relaxation  process  behind 
shock  waves.  The  thermodynamics  and  chemical  kinetics  included  in  the 
numerical  solutions  are  discussed.  Typical  results  obtained  from  the  normal 
shock  program  are  presented. 

Two  bow-shock  solutions  are  presented:  a  typical  high-altitude  re-entry 
condition  (i.  e.  M  =  22  at  200,  000  feet  altitude)  and  a  hypersonic  tunnel  test 
configuration.  The  rapid  expansion  of  the  flow  from  a  high-temperature 
reservoir  in  hypersonic  test  facilities  may  introduce  nonequilibrium  phenomena 
in  the  expansion  nocsle.  These  effects  will,  in  turn,  have  an  influence  on  the 
shock  layer  around  a  model  tn  the  expanded  flow.  A  finite-rate  expansion  solution, 
developed  at  Cornell  Aeronautical  Laboratory,  was  used  to  obtain  the 
expanded  flow  conditions.  These  nonequiUbrium  results  were  then  used  as  inputs 
to  the  bow  shock  solution.  For  both  cases,  the  calculation  was  carried  approximately 
65-  oft  the  centerline  and  the  computed  body  shapes  closely  resembled  a  circle. 

The  results  tor  the  re-entry  case  indicate  the  entire  shock  layer  to  be  in  a 
state  of  nonequilibrium  at  a  high  translational  temperature.  Results  for  both 
examples  are  presented  in  graphical  form  and  are  discussed  in  detail. 
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INTRODUCTION 


The  purpose  of  this  research  is  to  obtain  a  comprehensive  and  quantita¬ 
tive  understanding  of  the  infrared  radiation  from  the  heated  air  surrounding 
re-entry  bodies.  The  program  involves  the  detailed  spectral  measurements 
he  pertinent  radiating  species  and  an  analysis  of  the  aerodynamic  flow 
around  the  body.  This  analysis  takes  into  account  the  finite-rate 
chemistry  of  the  heated  air,  and  provides  the  concentration  profiles  of  the  gaseous 
constituents  and  the  local  thermodynamic  properties  in  the  flow  field. 

Three  earlier  progress  reports'-^  contain  details  concerning  the  sub¬ 
division  of  the  overall  project  into  its  respective  phases,  the  considerations 
which  have  determined  the  method  of  attack  on  each  phase,  and  a  description 
of  the  equipment  involved  on  all  phases.  Consequently,  these  items  will  be 
discussed  here  only  to  an  ea.ent  consistent  with  general  completene-ss. 

The  program  has  been  divided  into  two  areas.  The  first  involves  quanti¬ 
tative  measurements  on  the  system  of  radiating  species  of  air  at  elevated  tem¬ 
peratures;  a  machine  computation  of  the  state  of  the  air  in  the  flow  field 
surrounding  the  re-entry  bodies  constitutes  the  second  phPse.  The  con  b-'ration 
of  these  two  phases  will  allow  the  composite  radiation  picture  to  be  calculated 
tor  re-entry  systems.  A  portion  of  the  second  phase  involves  the  study  of  the 
chemical  kinetic  rates  necessary  to  couple  the  complex  chemistry  of  the 
heated  air  with  the  aerodynamic  flow  problem. 
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The  question  of  the  applicability  of  these  results  to  nonequilibrium  flow 
configurations  warrants  a  brief  discussion.  The  results  of  the  flow-field  analys 
presented  in  Section  H  of  this  report  show,  tor  example,  that  the  gases  in  the 
entire  region  between  the  bow  shock  and  the  blunt  body  may  be  in  a  nonequili¬ 
brium  distribution.  This  is  to  say  that  the  concentrations  of  the  various  species 
which  make  up  the  heated  air  are  not  those  which  would  be  obtained  from  a 
chemical-equilibrium  analysis  at  the  local  temperature.  This  is  due  to  the 
finite  rate  at  which  the  reactions  occur  in  the  flow.  However,  given  the  species 
concentrations  and  the  degree  of  excitation  of  the  molecules,  the  prediction 
of  the  spectral  radiation  can  still  be  reliably  made  if  the  transition  probabilities 
for  the  pertinent  transitions  of  the  optically-active  molecules  are  known. 

These  consideration,  gave  rise  to  the  basic  approach  made  upon  the  problem, 
which  renders  the  solution  applicable  to  the  nonequilibrium  problem. 

This  research  doe,  not  at  present  treat  the  case  of  nonequilibrium  radia¬ 
tion  arising  from  those  species  formed  within  the  shock  front  in  high-energy 
states.  A  general  experimental  approach  to  that  case  cannot  be  simply  made, 
since  a  complete  duplication  of  ambient  conditions  of  density,  composition,  and 

velocities  are  required.  That  problem  will  receive  special  attention  in  later 
phases  of  the  work. 

Section  1  of  this  report  deals  with  the  radiation  phase  of  the  program. 
Experimental  studies  were  performed  on  two  separate  problems  involving  radia¬ 
tion  from  heated  air.  One  was  the  continuum  radiation  associated  with  the 
recombination  of  NO  and  O  atoms  into  NO^.  This  continuum  extends  from 
about  4000  Angstroms  to  1.  5  microns  in  the  near  infrared.  Previously  reported 
measurements--’  were  confined  to  the  visible  region  of  the  spectrum,  and  the 


is 


2 


QM-  1373-A-4 


measured  rate  constants  were  obtained  at  room  temperatures.  The  experiments 
reported  herein  were  designed  to  measure  the  absolute  intensity  of  the  radiation 
in  the  near  infrared,  and  to  obtain  a  rate  constant  for  the  recombination  reaction 
which  applies  at  high  temperatures  (  3000°K). 

The  second  problem  was  the  measurement  of  the  absolute  intensity  of 
radiation  from  nitrogen  at  high  temperatures.  The  first  positive  band  system 
in  the  spectrum  of  nitrogen  extends  from  the  visible  through  the  near  infrared. 

In  this  wavelength  region,  it  is  the  strongest  of  the  radiating  species  in  heated 
air.  The  spectra  of  nitrogen  over  a  range  of  temperature  and  density  are 

presented.  The  method  of  data  reduction  are  also  described  and  the  results  are 
discus  sed. 

Section  II  discusses  the  machine  computation  programs  used  to  calculate 
the  state  of  the  gas  surrounding  re-entry  bodies  and  behind  normal  shock  waves. 

A  general  discussion  of  the  thermodynamics  and  chemistry  involved  in  the 
programs  is  presented.  The  finite-rate,  normal-shock  solution  has  been 
completely  programmed  for  an  IBM-704  computer,  and  has  a  capacity  for  20 
species,  40  chemical  reactions.  A  great  variety  of  test  cases  have  been  run. 
Results  are  presented  for  a  typical  air  calculation  using  a  s even- spec ie s ,  six- 
reaction  kinetic  model.  The  effect  of  coupling  the  degree  of  vibrational  excitation 
with  the  chemical  reactions  is  shown. 

The  inviscid  bow-shock  problem,  employing  the  inverse  method  of 
solution,  has  been  programmed  for  the  704  computer,  utilizing  the  large-capacity, 
high-speed  core  storage  unit.  At  the  present  time,  this  program  has  a  capacity 
for  1 0  species- and  20  chemical  reactions.  Several  test  cases  have  been 
completed  with  this  program,  for  both  a  diatomic  gas  and  air.  Results  for 
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a  typical  re-entry  vehicle  at  an  altitude  of  200.  000  feet  and  a  hypersonic  tunnel 

configuration  are  presented,  using  the  same  air  model  employed  in  the  normal 
shock  solution. 

It  IS  known  that  viscous  phenomena  may  effect  the  flow  field  around  a 
re-entry  body  in  a  rarefied  atmosphere.  However,  the  viscous  effect  is  not 
important  at  the  outer  portion  of  the  shock  layer  for  altitudes  of  approximately 
200,  000  feet  and  below.  The  ratio  of  boundary-layer  thickness  to  shock-layer 


thickness  at  the  stagnation 
this  upper  altitude.  Thus, 


region  of  hypersonic  vehicles  is  about  one-tenth 
the  inviscid  computation  can  be  used  with  some 


at 


degree  of  confidence  to  predict  the  chem i <! t -ptr  of  fv.ra  u  i  i 

p  cuieL  me  cnemistry  of  the  shock  layer  at  an  altitude 

of  200,  000  feet. 
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SECTION  I  -  RADIATION  MEASUREMENTS  PHASE 


The  „eed  lo.  .he  study  ol  .he  individual  components  in  the  determination 
Of  the  radiative  properties  of  a  multicomponent  gaseous  mixture  such  as  air  a. 
elevated  temperatures  has  been  stressed  in  Bte  previous  progress  reports  ^ 
Because  the  concentrations  of  the  radiating  species  each  vary  marhedly  and 
uniquely  With  Changes  in  the  temperature  and  density  of  the  gas,  a  prediction 
Of  the  radiation  from  the  composite  gas  under  arbitrary  conditions  requires  a 
hnovdedge  Of  the  detailed  radiative  behavior  of  each  component.  Therefore,  our 
P  P  as  been  to  assemble  the  necessary  apparatus  and  to  perform  detailed 
measurements  of  the  radiat.on  sp.Cra  in  order  to  enable  the  transition 
probabilities  for  the  component,  to  be  deduced. 

The  use  of  the  .hock  tube  a.  a  means  for  heating  gases  for  quantitative 

spectroscopic  studies  ha.  become  widespread*' ^  and  stem,  from  the  distinct 

advantage  that  gases  so  treated  are  brouglsto  thermodynamic  equilibrium  a. 

conditions  otherwise  unaffainable  in  ,h.  laboratory.  The  price  for  this  advan- 

.agc  lie,  in  ihe  oatremely  shor.  testing  time,  during  which  the  heated  gas  can 

he  Studied!  hence,  there  is  imposed  the  requirement  for  high-speed  spectre 
scopic  equipment. 

Complete  di.cti.sion.  of  criteria  lo  he  me.  by  both  shock  tubes  and 
apcc.ro.raph.  for  ,he  e.periment.  of  this  research  program  have  been  presented 

.n  Reference.  1  fhrou.h  3.  In  these  reports  the  details  of  the  completed  shock 

tube  and  the  tw^dve  -  channel  irfrared  spectrometc*  hp,  k 

-a  speciroinetei-  have  been  presented,  as 

well  a.  Uie  Hide  of  thia  combiimtion  in  soar.e  preliminar 

preiinninary  measurements  of  the 

fir*t  fxiailive  band  ayatem. 
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A. 


THE  NO-O  recombination  RADIATION  MEASUREMENTS 

I,,- 


^  The  Metory  of  th,  peoblen^  some  exploratory  measurements,  and  pro- 

^.mtnary  ata  were  presented  in  Reference  I,  the  last  Semi-Annual  Report. 

or  general  completeness  some  portions  of  that 

portions  of  that  report  will  be  included  in 

the  present  discussion. 

4ooo.K,Trrf 

K)  IS  that  of  the  radiative  recombination 

k 


NO  +  O 


NO^  +  hV 


reaction  has  been  ^-eported  by  Kaufman^.  It  was 

deduced  from  measurements  of  th^  • 

cinents  oi  the  intensity  of  the  a 

y  oi  the  radiation  from  discharge- 

generated  O  combining  with  NO  in  a  flo,  .  k 

ube  at  room  temperature  The 
reported  rate  is  =  I.  5  x  10^  cm^  mole'^  sec'^  Wth 

of  a  factor  of  five  Anoth  "" 

scuff  c  f  Ronfifn  and 

conf.rm  dUs  room-icmpcraiure  rsac.ion  rate  i„  an  accuracy  pf 

A,.p.-cciah,c  conccn.ra.ions  of  ihs  rs.canis  could  ho  .onorafed  in  fhe 

m-  u  a  moasuromon. 

-MO  consian.  could  ho  made  af  Mphor  iomporaiuros  and  pos.rhiy  wi.h 

.or  procsion.  I„  . . .  iomporaiuro  dopondonce  of  iho  raio  cons.anf 

might  be  deduced. 

2£i2£i:ii_Disc^^^  of  the  Experim^n^ 

HO,  formauon  has  boon  g.con  roconfiy  hy  Carom,  Owy„,  and  Mo.ho.iis  ^ 

-y  roporf  a  con.muum  whrch  hcg.ns  af  ahou.  dhOOA,  rcma.ns  falHyfiaf  ' 
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spectrum  out  to  1  micron,  although  no  intensity  measurements  were  made. 
Stewart  ^  has  published  a  profile  of  the  ultraviolet  end  of  this  continuum  in  arbi¬ 
trary  units.  The  first  quantitative  results  on  the  intensity  were  obtained  by 
Kaufman^  and  more  recently  by  Fontijn  and  Schiff^  who  measured  the  spec¬ 
trum  only  to  6200A.  The  sharp  cutoff  of  the  continuum  radiation  at  about  3600A 

from  the  radiative^  recombination  of  NO  and  O  at  room  temperature  has  been 

shown  by  Stewart  "  to  verify  Gaydon-s  analysis  of  the  energy  of  separation 

of  the  NO-O  unbound  system  from  the  ground  state  of  the  NO^  molecule.  In 

measurements  with  heated  gases  it  would  be  expected  that  the  additional  kinetic 

energy  of  the  reactant  species  would  cause  a  slight  shift  in  the  cutoff  wavelength 
into  the  ultraviolet. 

The  feasibility  of  making  a  quantitative  study  of  this  system  was  inves¬ 
tigated.  Since  these  experiments  were  performed  with  air.  «ear-maximum 
radiation  should  occur  at  the  greatest  value  of  the  NO-O  concentration  product. 
For  air  at  these  conditions,  the  NO^  concentration  is  small  (  ^  1  part  in  10“^  at 

3000»K),  and  is  maintained  at  its  equilibrium  values  chiefly  by  the  three-body 
reaction 


NO  +  O  +  M 


NO^  +  M 


The  rate  constant  for  this  reaction,  measured  at  room  temperature^®,  is 

^  °  mole  sec  .  The  value  of  k^  measured  at  room  tempera¬ 

ture  is  1.  5  X  lo"  cm^  mole-i  sec'^  With  these  two  values,  the  approximate 
mation  of  NO^  by  each  of  the  two  reactions  can  be  calculated.  Using 
the  equilibrium  values^  ^  for  the  concentrations  of  NO  and  O  in  high- to mpe nature 
air,  one  obtains,  at  .}000'’K 
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Radiative  recombination 


Tftree-body  recombination 


=  1.3x10 


-5 


(%) 


3/2 


moles  cm  ^  sec'^ 


d 

■HT 


[no,  ]  = 


1.  3 


moles  cm''*  sec 


where  the  approximate  densitv  j 

At  4000»K  each  of  the  production  re.^rnbouVCt 

cases  the  three-body  rate  i.  about  I0^imes  fa  ^ 

^-ation  The  three-body  reaction  Keep.  -tliative  recom- 

with  the  NO  and  Or-  ^  oncentration  in  equilibrium 

the  blackbody  value.  >-adiation  density  is  far  below 

riment, 

■A  Bausch  and  Lomb  lars«  Titters 

adapter’^  which  alio  •P®ctrograph  was  equipped  with  an 

Which  allows  four  photomultiplier-  to  K- 

of  the  spectroRraph  The  d-  Po**tmned  in  the  focal  plane 

An«,.ro„.,.  Thee  wa  ,  ""  *' 

. . . . 

discerned.  ultraviolet  .hould  have  been  readily 

.pact..,,.  m7in77'  ^  '--e 
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The  wavelength  distn.bution  eer,esp„„d=d  the  profile  that  had  been  reported 
However,  the  nae  of  this  ul.ravio!et-to-gree„  region  of  the  apec.run,  for 

tiuantitativ.  . . .  to  be  intpractica.  in  a  shock  tube  analysis. 

Th.s  i,  due  to  «.e  f.c  that  «  wavelength,  the  emission  from  the  Schumann- 

Rung.  system  of  tte  heated  O,  which  is  present  in  the  equilibrium  sample  of 
heated  air  in  the  tube  has  a  eisabl.  contribution.  This  was  determined  by  making 
etparate  runs  in  pure  oxygen  under  conditions  which  duplicated  the  oxygen 
content  in  the  heated  air  run.  The  intensities  representing  the  difference  between 
those  of  air  and  pure  oxygen  had  a.,  proper  distribution  as  reported  for  the  NO 
continuum.  It  was  apparent,  however,  that  the  errors  involved  in  such  a  sub-  " 
tractive  procedure  precluded  exact  measurements  and  also  obscured  any  effect 
such  a.  shift  Of  cutoff  wavelength.  Experimeata  in  pure  nitrogen  revealed  a 
totai  abaenc.  of  radiation.  m.r.by  eliminating  any  band  system  of  or  any  i„p„,i. 
.tee  a.  contributor,  to  the  mea.ured  radiation  in  th.  air  axperiments. 

Sev.ral  other  experimental  ech.me.  were  used  in  further  experiments 
•O  estabiieh  me  nature  of  the  epacrum  from  heated  air.  One  of  these  was 
the  use  of  a  Dumont  K-i430  phototube  and  a  Wratten  cutoff  filter  w.thou,  a 
dtsperaing  element,  Wim  this  „.,em  the  total  rad.aticn  between  0.  6  and 
1.  1  micron.  »a.  recorded,  and  the  re.ultan.  intensity  was  correlated  with  the 
reactant  concentrattons.  I.  b.can,.  apparent  that  to  analyse  the  radiation  quan- 

•..attvcly,  tb.  mtansity  have  to  be  obtained  sp.ctrallyi  i.  e.  .  as  a  function 

of  vvfilvc a.’ 

e  entire  spectruin  mapped  for  a  f^iven  set  of 

Utetmodynanue  condition..  accurate  value  of  the  rat,  constant  could  be  obt.mcd 
The  infrared  (>  1.  f 

cT.nd  Since  the  determination  of  the  rate  rnn«Ta«f-  a. 

a  e  constant  requires  the  iXicaBurement  of  the 

total  radiation,  it  was  felt  tVism  a-  .  . 

that  this  portion  of  tjie  spectrum  should  he  examined. 
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This  was  conveniently  done  by  making  use  of  the  multichannel  infrared 
spectrometer  which  had  already  been  used  in  the  early  worl  A  series  of 
experiments  showed  a  slow  drop  in  the  intensity  of  the  light  at  .  nger  wavelengths, 
which  between  1.  5  and  1.  8  microns  decreased  to  about  10%  of  its  value  at  1  micron 
It  was,  however,  impractical  to  proceed  with  the  overall  spectrum  measurement 
v/ith  this  instrument  because  of  its  low  detector  response  at  0.  6  microns,  and 
because  the  twelve  channels  could  not  be  made  to  subtend  more  than  0  3  microns 
at  any  given  setting.  It  was  decided  to  avoid  the  procedure  of  shifting  the 
instrument  across  the  0.  5  to  1.5  micron  range,  while  trying  to  repeat  the  shock 
tube  conditions  identically  for  at  least  three  runs.  The  system  to  be  described 
was  therefore  designed  and  used 


Instrumentation  for  Total  Spectrum  Measurement 

The  windows  on  both  sides  of  the  shock  tube  were  used.  One  side  had 
the  multichannel  infrared  spectrometer  deployed  in  exactly  the  same  manner  as 
shown  in  Figures  3  and  5  of  Reference  3  The  spectromete  r  was  used  in  first 
order,  and  !iie  twel'^o  rh'.nir-lK  the  v  Ta!.;-e  from  G,  9  In  ]  Z 

microns  Five  of  the  twelve  channels  were  used  in  the  c n pe  rim ents . 

Ihe  window  on  the  other  side  of  the  tube  was  viewed  by  li^e  Dumont  K- 
1130  infrared-sciisitivc  phuLornull  iplie  r  tubes,  each  of  whicli  was  equipped  with 
a  B.iuscli  and  Lomii  into  rLerriico  lilter  and  an  appropriate  W  ratten  <  utofl  filter, 
so  that  eaeh  tube  ia-spond*  d  t»  very  narrow'  v/a  ,elenot|i  band  Th,'  half- 
widliis  at  the  niters  ranged  'letwcen  iOdand  3  10  Angstroms.  I’lm  filte  r  -  photornu  )ti  ,d. 
combinations  v-erc  nlactd  so  that  eai  ii  ol  them  saw  o.lmost  the  sumo  volume  of 
ga-.i.  Eaeh  tube  suldcnded  a  viewing  cone  defined  by  <ui  apeiturc  of  1-inch  diameter 
at  a  distance  of  about  IM  mt'nes  from  ihe  shock  tube  wmdow.  llie  five  overUpping 
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cones  were  contained  in  an  aperture  of  2  inches  at  a  distance  of  12  inches. 

This  was  accomplished  by  means  of  a  five -sided  pyramidal  mirror  assembly 
which  reflects  the  total  beam  into  five  directions  in  a  radial  manner  (see  F'igui  e  i) 
The  filter-photomultiplier  assemblies  are  mounted  on  a  14-inch  diameter 
circle.  The  calibration  of  this  system  was  accomplished  in  the  same  manner 
as  that  of  the  multichannel  spectrometer,  namely  by  effectively  replacing  the 
image  of  the  shock  tube  window  with  that  of  a  calibrated  tungsten  ribbon- 
filament. 

The  five  wavelengths  chosen  for  study  are  given  in  the  following  table  in 
addition  to  those  covered  by  the  multichannel  spectromete** 

Instrument  Waveleng  th  (Microns) 

.  542 
.  598 

Photomultiplier  Ensemble  695 

798 

.  906 


.  906 
955 

Mulliclianncl  Infrared  Spectrometer  1.  005 

1  f)85 
1  190 

The  .  906  rnicvon  wavelength  was  recorded  by  both  systems  and  si.m  ^  e;! 
ae  a  cronscheck  on  the  overall  calibration  aifuracy  and  systerri  responsc- 

Each  shock  tube  rui*  conelitutcd  a  cov'crage  of  the  continuum  ai  Ib.c  nine 
wavelength*  •ho'a.n.  The  radiation  intersities  were  displayed  as  a  ]\inc»ic>n  o! 
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time  at  a  sweep  rate  of  20  microseconds/cm.  In  addition,  the  dual  beam 
oscilloscope  which  recorded  the  .906  micron  channel  using  the  two  systems  was 
also  coupled  to  another  oscilloscope  which  recorded  at  a  slower  rate  (200 
microseconds/cm)  and  served  to  monitor  the  long-time  radiative  behavior  of  the 

shocked  gas.  Reflected  shock  pressures  were  also  recorded  to  monitor  shock- 
tube  performance  for  each  run. 

The  general  features  of  the  radiation  records  included  a  rise  of  radiation 
as  the  reflected  shock  wave  swept  by  the  cone  of  gas  viewed  by  the  respective 
channels,  followed  by  a  reasonably  conatant  radiation  plateau.  Thi.  constant 
level  of  radiation  in  each  channel  indicated  the rmodynaniic  eqoiUbtium,  and 
constituted  the  datum  for  the  inten.ity  et  that  vavelength.  After  each  experiment, 
the  calibration  lamp  radiation  at  the  pre.cribed  current  value,  waa  pa.sed 
through  chopper,  and  pre.eMed  to  the  detectora.  These  intensitie.  were  also 
displayed  and  photographed.  Thua,  a  compariaon  of  Ute  signal,  from  .hock 

tube  and  lamp  yielded  abaolut.  value,  for  the  inten.ity  of  the  radiation  from  the 

shock-heated  gas. 

The  results  of  the  first  experiments  with  this  apparatus  were  presented 
and  discussed  in  Reference  1.  The  experiments  were  m  progress  at  that  time 

e-ind  the  results  were  presented  as  preliminary.  They  will  be  summarized  at  this 

point. 

Several  spectra  were  obtained  which  showed  in  general  the  form  of  the 
recombination  radiation  spectrum.  The  intensity  m  the  nrne  channels  wa.s 
reported  in  watts  cm"^  ste  radian'^  micron'  ^ 

The  probable  error  m  the  data  points  due  to  the  readability  of  the  records 
uas  about  1%.  In  addition,  the  largest  difference  in  the  radiation  intensity  at 
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.  906  micron  as  recorded  by  the  two  separate  detecting  systems  was  7%,  while 
for  most  runs  the  agreement  wa«  within  3%.  This  difference  was  an  indication 
of  the  precision  of  the  calibration  procedure,  since  both  tungsten  lamps  were 
calibrated  against  the  same  NBS  standard  lamp.  The  largest  source  of  error 
in  the  results  arose  from  the  problem  of  determining  the  best  time  (measured 
from  shock  wave  reflection)  at  which  to  measure  the  radiation  intensity  level. 
Variations  between  3%  and  10%  were  sometimes  recorded  over  the  testing  time 
interval,  and  hence  gave  rise  to  an  uncertainty  in  the  absolute  spectral  intensity. 
The  general  procedure  was  to  record  the  level  of  intensity  at  about  120  to  140 
microseconds  from  the  time  of  shock  wave  reflection.  The  reflected  shock  wave 
required  some  30  microseconds  to  move  through  the  gas  in  the  field  of  view,  and 
an  additional  100  microseconds  was  ample  time  for  complete  equilibration  in  the 
tube.  The  total  testing  time,  measured  by  both  pressure  and  radiation  channels, 
ranged  between  150  and  300  microseconds. 

There  were  several  indications  that  the  origin  of  the  radiation  was  the  NO 

id 

rer  mbination.  Primarily,  over  a  broad  range  of  conditions  the  radiation 
intensity  could  be  normalized  by  the  product  of  the  NO  and  O  concentrations. 

In  addition,  check  runs  were  made  in  pure  and  O^.  For  nitrogen  at  T=4240“K 
at  4.  25  times  normal  density,  no  radiation  was  recorded  in  any  channel.  This 
indicated  that  the  source  of  the  radiation  must  be  linked  with  the  presence  of 
oxygen,  and  further,  that  no  radiation  shock-tube  impurities  such  as  CN  are  present 
in  these  experiments.  The  air  used  in  the  present  studies  was  procured  as  "dry 
air"  from  the  Matheson  Company.  One  check  on  the  results  at  high  tcmperaiures 
was  mad*  by  performing  a  run  using  an  air  mixed  from  reagent  grade  and  O 
The  results  were  essentially  the  same,  and  indicated  no  noticeable  impurities  in 
the  dry-air  container. 
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Experiments  in  which  pure  oxygen  was  brought  to  the  same  conditions  as 
those  of  the  high-temperature  air  runs  yielded  only  a  very  small  amount  of 
radiation.  It  was  almost  immeasurable,  due  to  the  low  signal  -  to  -  noise  ratio. 

It  should  be  noted  that  experiments  in  air  and  various  mixtures  have  been 

included  in  the  schedule  of  studies  to  be  performed  directly  after  the  measure¬ 
ments  on  the  first-positive  band  system. 

Preliminary  estimates  of  the  rate  constant  were  made  from  the  data, 
and  revealed  that  the  rate  at  about  3500‘K  rvas  some  20  times  lower  than  the 
room-temperature  rate  reported  by  Kaufman^  It  was  noticed  in  the  spectra, 
however,  that  at  0.  6  microns  a  large  spike  of  radiation  was  recorded,  for  experi¬ 
ments  in  excess  of  4000°K.  Because  of  the  ’  kelihood  of  radiation  from  sodium 
at  this  wavelength  it  was  decided  to  place  the  interference  filters  (which  determine 
the  bandpass  of  the  photomultiplier  detector  channels)  at  wavelengths  which 
would  be  free  of  radiation  from  such  unavoidable  impurities  as  sodium,  calcium, 


and  potassium. 

The  final  placement  of  the  filters  was  determined  after  first  photographing 
as  much  of  the  continuum  spectrum  as  possible.  A  medium  quarts  spectrograph 
(f  lO)  was  used  together  with  Kodak  1C3-F  and  I-N  spectrographic  plates  to 
obtain  the  spectra  shown  in  Figure  2.  The  first  spectrum  is  Uial  for  air  at 
about  4500 ”K  and  J  times  normal  density,  in  which  the  plate  was  exposed  to  the 
total  shock  Mbe  radiation  phenomenon.  Tliis  includes  not  only  the  radiation  from 
Uie  gas  during  the  "clean”  testing  time,  but  also  the  subsequent  radiation,  largely 
from  impuriUes  and  from  the  interaction  of  the  hyd  rogen -ci  r  i  ve  r  gas  with  the 
heated  tost  air.  The  prominent  OH  band  at  3064  Angstroms  is  clearly  seen 
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The  use  of  a  capping  shutter  which  limits  the  plate  exposure  only  to 
the  radiation  from  the  gas  after  processing  by  the  reflected  shock  permitted  the 
second  spectrum  of  Figure  2  to  be  obtained.  The  radiation  from  the  Schumann- 
Runge  band  system  of  oxygen  can  be  clearly  seen  on  the  plate.  The  radiation 
of  interest  to  the  present  problexxi  extends  from  about  4000  Angstroms  into  the 
infrared.  The  prominent  lines  in  this  region  are  those  of  potassium,  calcium, 
and  sodium.  The  unexposed  region  between  5000  and  5500  Angstroms  results 
from  a  decrease  in  the  sensitivity  of  the  emulsion.  The  type  F  plates  were 
used  to  examine  this  region  in  detail.  The  last  spectrum,  together  with  the 
other  results,  show  that  this  radiation  does  not  appear  in  the  spectra  of  either 
pure  oxygen  or  nitrogen. 

The  interference  filters  were  placed  in  the  cleanest  portions  of  this 
continuum,  and  were  now  at  the  wavelengths:  4750,  5100,  6700,  8100,  and 
9060A.  With  this  configuration,  another  series  of  experiments  were  performed 
with  air  at  temperatures  between  2800  and  4500*K,  and  densities  between  2 
and  20  times  normal. 

Results 

The  data  from  a  series  of  experiments  near  3700 °K  are  shown  in  Figure 
3.  In  this  graph  the  measured  intensities  at  the  various  wavelengths  are 
plotted  against  wavelength.  The  intensity  was  normalized  by  the  product  of  ibe  NO 

and  O  concentrations; 

where  Z  is  the  compressibil  ity  f  actor,  Y*  fhe  air  density,  and  the  X's  arc 
mole  fractions,  obtained  from  Reference  11.  Tpe  smooth  curve  of  the  graph 
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between  0.  4  and  0.  62  microns  was  taken  from  the  detailed  spectrum  by 
Fontijn  and  Schiff  .  which  was  reported  during  the  course  of  the  present  work. 
The  spectrum  was  fitted  to  the  data  at  0.  51  microns. 

As  can  be  seen,  the  agreement  with  the  intensity  distribution  verifies 
the  recombination  radiation  as  the  chief  contributor  to  the  recorded  spectrum. 
The  experiments  covered  a  range  of  3  in  the  factor  ^  ^ 

the  intensities  are  proportional  to  it  to  within  20%. 

The  rate  constant  derived  from  these  data  at  3750°K  for  the  reaction 


NO  +  O 


NO^  +  hv* 


is  equal  to  4.  4  X  10^  cm^  mole'*  sec'*  +25%.  This  is  about  20  times  smaller 
than  the  room-tempe rature  value. 

The  reaction  rate  measured  at  is  lower  than  that  reported  for 

300”K  by  a  factor  of  about  (3750/300)' ^  which  is  reasonable  for  the  temperature 
variation  of  recombination  reactions.  However,  between  3000  and  3750°K  the 
values  show  a  T'^'  ^  dependence.  This  variation  cannot  reasonably  be  attributed 
to  a  change  in  the  reaction  rate  constant,  but  rather,  indicates  a  contribution  of 
radiation  from  another  source  at  these  lower  temperatures. 

In  contrast  with  the  room-temperature  work,  for  air  in  equilibrium  at 
these  conditions,  radiation  can  arise  both  from  the  NO  -  O  recombination  and 
from  excited  states  of  equilibrated  NO^  molecules.  The  density  dependence  of 
the  radiation  from  these  two  sources  is  the  same.  However,  the  equilibrium 
population  of  the  NO,  falls  rapidly  to  a  negligible  value  by  3  750'’K*'*.  Further, 
subtraction  of  the  NO-O  recombination  intensity  at  the  lower  temperatures  (based 
on  a  T  dependence)  leaves  a  residue  which  ia  consistent  with  the  radiation 
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from  NO^  based  an  estimates  made  from  the  absorption  coefficient  of  the 
NO^  molecule.  A  recent  paper  by  Broida  and  Schiff  and  Sugden  shows  a 
structured  spectrum  at  the  cutoff  of  the  room  temperature  NO-O  recombination 
continuum  which  suggest*  the  presence  of  excited  NO^  molecules  in  their 
experiments.  The  precision  of  our  data  doe*  not  warrant  further  analysis  of 
the  NO^  contribution. 

One  experiment  was  performed  at  2800*K.  At  these  temperatures  the 
radiation  was  seen  to  rise  gradually  in  time  to  an  equilibrium  value  (for  higher 
temperatures  the  rate  of  rise  was  due  only  to  the  transit  time  of  the  reflected 
shock  wave  across  the  viewing  port).  The  rate  was  strongly  suggestive  of  the 
relaxation  time  to  chemical  equilibrium.  A  check  was  made  on  the  kinetics 
of  this  shock  tube  run  by  using  the  normal  shock  program  discussed  in  Section  I! 
of  this  report.  These  results  gave  good  agreement  between  the  measured  rise 
of  the  radiation  and  the  formation  of  NO  molecules  which  is  rate-^determining 
for  NOj  formation.  This  experiment,  coupled  with  the  absence  of  radiation 
in  pure  or  experiments  serve  to  confirm  the  source  of  the  radiation 
and  the  validity  of  the  results. 

B.  THE  NITROGEN  INFRARED  MEASUREMENTS 

One  of  the  principal  radiating  species  in  heated  air  is  the  nitrogen  molcc  ult 
whose  first  positive  band  system  extends  from  the  visible  to  about  two  mirrams 
in  the  infrared.  The  only  reported  transition  probabilities  for  this  system 
were  deduced  from  a  few  broad-band  measurements  made  in  heated  air  .  Ih<' 
need  for  a  detailed  study  of  this  system  was  dear,  and  thus  it  became  a  primary 
aim  of  the  present  research  program. 
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Description  of  the  Experiment 


The  apparatus  for  these  experiments  consists  basically  of  the  shock  tube, 
used  to  process  the  gas»  and  the  twelve -channel  infrared  spectrometer,  with 
^hich  the  absolute  intensity  of  the  spectrum  of  the  gas  was  measured.  The 
criteria  leading  to  the  design  of  these  components  together  with  the  details  of 
construction  and  calibration  have  been  reported  in  References  2  and  3.  Briefly, 
the  shock  tube  has  a  teat  section  30  feet  long,  of  3-inch  inside  diameter.  It 
18  equipped  with  contoured,  sapphire  windows  near  the  reflecting  wall.  It  is 
designed  for  operation  at  high  driver-gas  pressures  (30,  000  pei)  and  at  good 
purity  levels  (total  rate  of  rise  at  vacuum  ^  2  microns/hr.  ).  Combustion  drive 
was  not  used  because  of  the  inherent  difficulties  of  performance  and  cleanliness 
with  combustion-driven  shock  tubes. 

The  twelve- channel  infrared  spectrometer  was  designed  to  span  the 
range  from  0.  7  to  2.  4  microns.  It  i.  an  f/3.  5  Ebert-type  grating  instrument 
and  the  twelve  channels  subtend  about  0.  3  microns  in  first  order,  with  a  bandpass 
of  0.  01  microns  per  channel,  each  of  which  has  a  lead  selenide  detector  with 
about  10  microsecond  response  time. 

The  use  of  this  system  and  some  preliminary  measurement,  of  the  nitrogen 
system  has  been  reported  in  Reference  3.  Since  that  time  the  system  has  been 
modified  and  refined,  particularly  in  regard  to  the  radiation  calibration  tech¬ 
nique.  and  to  actual  shock  tube  operation.  These  will  be  described  ir  the  next  . 
subsection. 

Experimental  Procedure 

The  shock  tube  was  thoroughly  scrubbed  with  alcohol  and  hard-pumped 
between  runs.  Prepurified  nitrogen  from  the  Matheeon  Company  was  used  for  the 
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experiments.  The  following  mass-spectrometer  analysis  of  the  nitrogen  was 
obtained  from  the  Matheson  Company: 


^2 

H^O 

Hydrocarbons 

CO^ 

CO 

H., 


5  ppm 
3  ppm 

less  than  0.  5  ppm 
non-detectable 
non-detectable 
non-detectable 


The  gas  was  used  with  no  further  processing.  Initial  test  gas  pressures 
were  read  on  a  manifold  of  Wallace  and  Tiernan  absolute  pressure  gauges,  which 
were  calibrated  and  cross-checked  with  each  other.  The  speed  of  the  incident 
shock  wave  was  measured  with  thin-film  ionization  gauges  at  four  intervals  in 
the  last  six  feet  of  the  diock  tube.  The  intervals  are  18  inches  apart  and  the 

elapsed  time  was  measured  with  a  digitalized  output  unit  to  the  nearest  micro¬ 
second. 


In  order  to  traverse  a  spectrum  from  0.  7  to  1.2  microns  in  the  second 
order  {which  was  done  in  these  experiments),  the  conditions  of  the  gas  must  be 
reproduced  exactly  for  five  to  six  shock  tube  runs.  Using  only  pressure  breaks 
on  the  diaphragms  at  the  2500  psi  level  makes  this  almost  impossible,  due  to  the 
unavoidable  variations  in  the  diaphragms.  This  problem  was  overcome  through 
the  use  of  a  double -diaphragm  technique.  The  double  diaphragm  unit  consists 
of  two  diaphragms  separated  by  a  3-inch  spacer.  During  the  loading  process, 
this  small  chamber  is  automatically  naaintained  at  one-half  the  driver-gas 
pressure.  For  desired  operation  at  2500  psi,  the  chamber  contains  1250  psi. 

Since  the  driven  section  is  at  a  comparative  vacuum  (10-200  mm  Hg),  the  diaphrag 
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support  a  pressure  differential  of  1250  psi.  Upon  venting  the  chamber,  the  full 
2500  psi  bears  first  upon  the  diaphragm  closest  to  the  driver,  and,  when  it 
ruptures,  upon  the  second.  Thus,  the  only  requirement  upon  the  diaphragms 
IS  that  they  support  1250  psi,  and  break  at  less  than  2500  psi.  This  is  easily 
accomplished,  since  normal  variations  in  the  rupture  pressure  of  the  diaphragms 
are  less  than  +  200  psi.  With  this  system  a  series  of  runs  can  be  made,  in 
which  the  wave  speed  in  any  of  the  timing  intervals  is  duplicated  to  within  1 
microsecond  out  of  150.  This  reproducibility  is  within  the  error  in  the  determina 
tion  of  the  density  and  temperature  behind  the  reflected  shock.  Thus  the  data 
require  no  corrective  normalization  between  successive  experiments  in  a  given 
series  at  constant  reflected- shock  condition- 

The  absolute  intensity  calibration  of  the  spectrometer  was  carried 
out  in  the  following  manner.  During  the  shock  tube  run,  the  entrance  slit  (1  x 
12  mm)  of  the  spectrometer  is  focused  on  the  sapphire  window  of  the  shock 
tube.  After  the  run,  it  is  focused  onto  the  tungsten  ribbon  filament  of  a 
standard  lamp.  This  is  done  without  changing  the  optical  geometry.  A  fast 
chopper  is  employed  to  produce  a  square  wave  output  which  is  recorded  by  each 
channel.  Thus  the  primary  data  consist  of  two  records  which  are  the  ratio  of 
the  radiation  from  the  gas  and  from  the  tungsten  lamp. 

The  voltage  output  of  a  typical  channel  is  shown  in  the  following  sketch: 


MICROSECONDS 
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The  sharp  rise  of  the  signal  corresponds  to  the  transit  time  of  the 
reflected  shock  wave  through  the  subtended  volume  of  the  test  gas.  A  rather 
linear  rise  in  signal  follows,  and  corresponds  to  a  radiation  increase  which  is 
felt  to  be  due  to  a  small  rise  in  the  temperature  of  the  gas  in  the  tube.  The 
radiation  is  highly  temperature  dependent,  and  Che  temperature  rise  needed  for 
this  increase  is  too  small  to  be  correlated  with  the  reflected  shock  pressure, 
output  over  the  total  interval  remained  constant.  The  datum  for  each 
channel  was  taken  to  be  the  intercept  of  the  linearly  extrapolated  slope  at 

'  ^-Pe-ture  .ould  be  that  calculated  by  simple  shock  tube  theory. 

g  standard  lamps  were  obtained  by  comparison  with  a  National 

Bureau  of  Standards  Tamp  whose  output  was  measured  and  certified  at  NBS. 

The  comparison  was  made  over  a  broad  range  of  wavelengths  from  0.  5  to  2.  4 

crons  by  alternately  passing  the  radiation  from  the  NBS  lamp  and  the  working 

standard  lamps  into  a  Gaertner  monochromator  TKo 

ocnromator.  The  current  in  the  NBS  lamp 

was  provided  by  a  stabilized  power  suonlv  whir-K 

power  supply  which  was  operated  from  an  0.  01%  line 

regulator.  The  current  was  set  at  the  i  -  , 

et  at  the  NBS  value  with  an  0.  2  5%  Weston  370 

an  0. 
amp  .n 

.Kp  cai..„,i„p.  Both  lead  .o, ph.de  and  .nrnaned-.en.l.ltee  pho,o„.,l..p,.e n  ...hea 

we.e  used  in  the  calibration  to  maintain  the  highest  signal-to-noise  ratios  and  to 

minimize  the  reading  error  ir,  tv,... 

C  ei.or  ,he  co„p„..o„  p. 

Overall-system  linearity  was  checkerl  K, 

y  checkt,d  by  comparing  the  readings  for  the  NBS 

lamp  ladiance  at  25  and  30  am,«..res.  since  the  NBS  lamp  was  calibrated  at  both 
The  effect  of  aging  of  the  Ump.  was  studied,  and  the  determina- 
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tion  of  the  usable  working  areas  of  the  filaments  was  made  by  measuring  the 
radiance  gradients  along  the  filaments. 

The  accuracy  of  the  overall  calibration  of  the  working  standard  is  felt 
to  be  better  than  5%.  The  largest  error  in  the  absolute  values  used  in  the  experi¬ 
ments  is  that  reported  for  the  NBS  lamp  itself^^. 

The  spectrometer  was  checked  out  to  insure  that  no  radiation  from 
unwanted  orders  or  wavelengths  were  overlapping  the  detectors  at  any  given 
setting.  Use  was  made  of  silicon  windows,  water  cells,  Kodak  infrared  filters 
and  several  organic  liquid  cells  in  these  checkouts  and  in  the  runs  where  necessary 
In  general,  only  the  radiation  subtended  by  tlie  detector  array  was  permitted  to 
enter  the  spectrometer.  The  intensity  calibr-.tions  were  always  made  with  the 
order- sorting  filters  in  the  system. 

The  linearity  of  the  detecting  system  was  also  checked  by  separately 
verifying  the  linearity  of  the  electronic  system  of  amplifiers  and  by  an  inverse- 
square  experiment,  where  the  response  of  the  detectors  was  measured  as  a 
function  of  the  distance  from  a  point  source.  No  deviation  from  linearity  was 

noted  for  values  of  detector  output  which  greatly  exceeded  those  obtained  in  the 
experiments. 

Results 


An  exploratory  series  of  runs  using  the  multichannel  spectrometer  in  first 
order  revealed  that  the  prominent  feature,  of  the  spectrum  lay  between  the  low 
wavelength  limit  of  0.  7  microns  and  about  1.  5  microns.  This  region  of  the 
spectrum  was  then  carefully  .surveyed  in  second  order,  by  making  a  series  of 
run,  at  conditions  of  constant  temperature  and  den.ity.  The  resultant  spectrum 
is  shown  in  Figure  4,  where  Uie  radiation  i.  plotted  in  watts  cm'^  steradian*  ^ 
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.  The  twelve  data  points  from  a  single  run  are  shown  as  separate 

symbols,  and  the  precision  in  the  data  can  he  determined  hy  examining  the  spread 

m  mtensttT  a.  a  given  wavelength  for  a  number  of  overlapping  runs.  Also  plotted 

g^^e  abscissa  are  the  pertinent  band-head  positions  according  to  Dieke  and 

Heath  ,  with  the  heights  of  the  markers  scaled  to  the  values  of  the  overlap 

integrals  after  darmain  and  Nicholls>«  The  vibrational  band  sequences  can  be 

clearly  seen,  and  indicate  the  advantages  of  the  multichannel  approach,  in  that 

an  radiation  corresponding  to  a  given  symbol  on  the  plot  was  recorded  from  the 
same  sample  of  gas. 

The  radiation  from  the  standard  lamp  is  given  in  units  of  watts  cm'^ 
ste radian  micron-'.  Conversion  to  the  units  plotted  is  accomplished  by 

dividing  by  the  path  length  through  the  radiatin,, 

®  S  g^s,  equal  to  the  inside  diameter 

of  the  shock  tube  (7.  60  cm  was  usedl 

he  sapphire  window-transmission  correc¬ 
tion  was  taken  to  be  0  9  and  i 

.  9,  nd  finally,  the  data  were  corrected  by  a  factor  of 

1.  4  to  account  for  the  effect  of  reflections  inside  the  shock  tube.  This  factor 

was  Obtained  by  experiments  described  in  the  Appendix.  Thus,  the  intensity 

plotted  in  Figure  4  was  directly  obtained  from  the  data  by  the  relation 


1  {- 


Watts 


cm  -  ster-mit  ron 


-)  = 


T/.b)(.  9)(1.4)- 


(- 


watts 


cm  -  s  te  r  -  mic  ron 


'n  Figure  5  are  presented  several  spectra  spanning  one  o,  the  vibrational 

band  sequences  (V  =  1,  ....c,  of  this  first  positive  system  _>.  A^S  , 

between  0.  78  and  0.  91  micronc;  ti  „ 

■crons.  These  spectra  were  obtained  by  varying  the 

density  and  temperature  of  (hn 

gen  in  experiments  with  a  constant  wavelength 

setting  of  the  spectrometer.  The  essentit,!  ,  r 

essential  features  of  the  spectrum  are  the  same, 
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and  the  variation  in  intensity  results  from  the  change  in  the  population  of  the 
upper  energy  level  involved  in  the  transitions.  This  effect  is  more  clearly 
seen  in  Figure  6,  which  gives  the  temperature  dependence  of  the  radiation.  Con¬ 
sidering  the  radiation  to  arise  from  an  upper  level  of  energy  E  ,  the  intensity 
is  proportional  with  the  population  in  this  level; 


I  =  c 


-E/kT 


where  n  is  the  number  of  molecules  per  cm^,  L,  the  Loschmidt  number, 

and  c  IS  a  constant.  is  the  internal  partition  function  for  nitrogen.  It  is 

plotted  in  Figure  7,  and  was  obtained  from  Reference  11.  Dividing  and  taking 
the  logarithm,  one  obtains 


0 

(-/Lj 


(4-) 


Thus,  a  plot  of  the  left  side  vs.  1/T  should  yield  a  straight  line,  whose 
slope  IS  -  E/k,  and  whose  intercept  permits  the  constant  c  to  be  evaluated. 
Eight  temperatures  between  6000  and  7200“K  are  plotted  for  each  of  two  wave¬ 
lengths,  chosen  at  each  edge  of  the  recorded  band  sequence.  It  can  be  seen 
that  a  good  straight  line  can  be  passed  through  the  points.  The  slope  of  the  line 
at  0.  887  microns  which  contains  the  head  of  the  (1,  0)  band  corresponds  to 
E  =  58800  cm  .  This  is  in  very  good  agreement  with  the  best  reported  value 
of  61700  cm  ,  and  substantiates  the  radiation  origin  to  be  N^.  Keck  et  al*^ 
have  reported  the  f  number  for  this  band  system  in  nitrogen  to  be  0.025. 

The  results  of  the  present  work  are  a  factor  of  five  below  their  value.  This 
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discrepancy  is  sigpificai^t  and  is  illustrated  in  Figure  8.  Because  their  results 
were  deduced  from  radiation  data  taken  in  shcck-heated  air,  the  discrepancy 
may  ba  due  to  a  prominent  radiator  present  in  air  but  not  in  nitrogen.  This 
point  was  checked  by  repeating  the  0.  78  and  0.  91  micron  experiment  with  air 
at  T  =  5400-K  and  a  density  of  0.  ?  times  normal.  The  resultant  radiation  was 
normalized  to  the  amount  of  molecular  nitrogen  present  in  the  air  at  these 
conditions.  At  0.  8  microns  the  radiation  intensity  was  equal  to  that  predicted 
by  the  present  pure -nitrogen  work,  indicating  the  absence  of  other  radiation 
species  in  this  wavelength  region,  at  least  for  temperatures  below  6000°K. 

Thus  it  is  felt  that  the  lower  intensity  values  for  the  nitrogen  spectrum  measured 
in  these  studies  are  valid. 

The  spectral  resolution  of  the  data  allows  the  transition  probabilities 
to  be  determined  with  good  precision.  Because  this  analysis  is  presently 
being  made,  it  could  not  be  included  in  this  report,  but  will  be  submitted  for 
publication  shortly.  Based  on  the  calculations  of  Keck,  et  ai,  the  f  number 
should  be  about  0.  005,  one -fifth  of  their  reported  value. 


SUMMARY  AND  CONCLUDING  REMARKS  -  RADIATION  PHASE 
The  reaction  rate  constant  for  the  NO -O  radiative  recombination 


NO  +  O 

lias  been  measured  to  be  4.  4  x 
3750°K,  This  rate  is  lO  times 


k 

10"  cm  mole  sec  +Z5%  at  temperatures  of 
smaller  than  that  reported  for  this  reaction  at 


room  temperature. 
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« 

2.  At  lower  temperatures  (2800-3100»K)  additional  radiation  was  recorded 
and  is  attributed  to  the  bound  states  of  thermally  equilibrated  NO 

2 

3.  The  detailed  and  intensity-calibrated  spectrum  of  the  first  positive  band 
system  of  nitrogen  has  been  obtained  between  0.  74  and  1 .  2  microns  for  a  range 

of  temperatures  between  5800  and  7200»K.  The  observed  temperature  dependence 
of  the  radiation  is  in  good  agreement  with  equilibrium- radiation  theory  and  clearly 
identifies  nitrogen  as  the  radiating  gas. 

4.  The  intensity  of  the  radiation  from  pure  nitrogen  is  lower  than  previously 

9 

reported  values  by  a  factor  of  five. 

The  scope  of  the  continuing  research  in  this  area  includes  the  experimental 
determination  of  the  other  contributors  to  the  radiation  from  heated  air.  These 
are  presently  felt  to  be  the  continuum  radiation  arising  from  the  interaction 
of  electrons  with  the  atoms  and  molecules  of  the  heated  air.  It  is  planned  to  study 
the  radiation  for  a  broad  range  of  -O^  ratios,  from  which  the  effectiveness 
of  the  various  scatterers  may  be  deduced. 
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FLOW  FIELD  ANALYSIS  PHASE 


SECTION  IT  -  RE-EnTRY  BODY 


The  solution  of  the  inviscid  nonequilibrium  blunt-body  problem  (which  has 

.ee„  ae.e,e„ce3  ,,  Z  a„.  3,  .3  nee^eO  .He 

.story  of  rad.at.o„  from  the  ho.  gaa  surrounding  a  re-entry  vehicle.  The  present 

.bed)  and  employs  coupled  vibration-dissociation  reactions  to  describe 

.be  thermoohemioa,  state  of  the  gas  behind  the  bow  shoch  wave.  This  genera, 

bow  shooh  problem  has  been  programmed  for  an  IBM-,04  computer  1„  addition, 

s  similar  finito-rate  solution  has  been  programmed  for  normal  shock  waves 

.0  assess  the  variation  of  rate  constants  on  the  over-al,  relaxation  none.  .  general 

escription  of  both  programs  and  pertinent  results  obtained  .0  date  will  be 
discussed. 

t.  .a  known  that  the  flow  field  around  a  re-entry  vehicle  in  a  rarefied 

atmosphere  may  be  influenced  by  viscous  effects.  Cheng^»  defines  a  flow  regime 

w  em  boundary  layer  or  vor.ici.y  interaction  theories  may  apply  as  O0)<eK^<^ 

wheie  IS  the  inverse  density  ratio  over  the  shock  c  h 

me  snock.  K  can  be  approximated 
by 


A 


eO 


-  -e  mean  free  path  upstream  o,  f,  . 

re-entry  vehicle  (i.  e.  one  foot  nose  diameter),  at  an  altifcli 
of  200,  000  feet,  EK^-^soo 

■  this  parameter  is  about  50,  or  apprea.  n  ,. 

a.  ZhO,  000  feet.  Scala's^-  similar  results  show  the  ratio  of  boundary  layer 
tbickness  io  Shock  layer  thickness  in  the  stagnation  region  of  hypersonic  vehi,  irs 
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is  about  1/  10  at  200,  OOO  feet,  and  increases  to  i/2  at  250,  000  feet  altitude.  The 
inviscid  solution  can  be  used  with  some  degree  of  confidence  to  predict  the  chemistry 
of  the  shock  layer  at  200,  000  feet. 

A.  NORMAL  SHOCK  SOLUTION 

General  Discussion 

The  normal  shock  program  is  similar  to  the  bow  shock  problem  in  the 
numerical  aspects  of  the  solution,  and  is  identical  with  respect  to  the  thermodynamics 
and  chemistry.  Since  it  is  a  much  less  complex  program  than  the  bow  shock 
solution,  the  general  characteristics  of  relaxation  phenomena  behind  shock  waves 
can  most  easily  be  discussed  in  terms  of  the  normal-shock  solution.  Consider 
the  sketch: 


■SHOCK. 


- 1 1 1  n  I  n  1 1 1 1 - 

- 

o 

In  shock  coordinates,  the  fiow  is  advancing  into  the  shock  wave  at  velocity  and 

leaving  at  velocity  U  .  When  the  air  passes  through  the  shock  wave,  it  ie 
compressed  and  most  of  the  kinetic  energy  available  in  the  free  stream  is  con¬ 
verted  to  thermal  energy.  The  active  (i.  e.  translation,  rotation)  contributions 
to  the  internal  energy  of  the  gas  axe  very  quickly  adjusted  to  th»  equilibrium  values 
in  the  time  of  a  few  molecular  collisioae.  How«v»r,  th*  »onactive  degrees  of 
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freedom  4i.  e.  vibration,  dissociation  and  ionization)  adjust  much  more  slowly. 

Duff  has  investigated  the  chemical  dis sociation- recombination  process  behind 
normal  shocks  in  air  coirsidering  sdbrational  equilibrium  of  the  diatomic  species. 
The  addition  of  vibyatienally  nonequilib rated  diatomic  species  and  the  coupling  of 

this  vibrational  relaxation  process  with  the  chemistry  has  been  studied  by 
Teare. 

The  mathematical  description  of  this  relaxation  process  used  in  the 
cal  program  is  cbtained  from  the  solution  of  a  set  of  simultaneous  differential 
equations  to  obtain  the  density  (  p  )  pressure  (  .jjlV  ),  flow  velocity  {U],  con¬ 
centration  of  thej^^  species  (  V.  .  vibrational  energy  of  the  J  species 
(  €j  )  at  any  distance  (  p  behind  the  shock  wave,  (Detailed  analyses  of  the 
normal  and  bow  shock  problems  can  be  found  in  References  2  and  3  and  will 
not  be  presented  here  )  The  conditions  immediately  behind  the  shock  wave  can 
be  determined  in  a  number  of  ways,  depending  on  the  thermochemical  descrip¬ 
tion  of  the  problem  In  general,  each  species  (  MJ  )  is  described  thermodyna¬ 
mically  with  translation,  rotation,  vibration  (described  by  harmonic  oscillator 
relations  for  diatomic  species),  and  excited  electronic  states  contributing  to 
the  internal  energy.  The  translational,  rotational  and  electronic  states  are 
assumed  to  equilibrate  immediately  belund  the  shock  and  detm-mine  the  initial 
nonoquilibrium  temperature.  In  flow  regimes  where  the  vibrational  relaxation 
time  for  the  diatomic  species  is  much  less  than  tnc  chemical  or  ionization 
1C. taxation  times,  the  oiatomic  species  can  be  assumed  to  be  equililiratcd  v  ibrat  iona  1 1- 
with  the  three  preceding  degrees  ol  freedom.  However,  for  low-density  regimes, 
vibrational  relaxation  times  become  long,  and  the  diatomic  species  will  be  in 
a  state  of  vibrational  noncquilibr lum  immediately  behind  the  shock.  This  is 
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(Jesct-ibed  tikrougja  a  system  of  harmonic  that  relax  to  equilibrium 

conditions  while  maintaining  a  Boltafnann  diatriJjutiori  about  a  vibrational  tempera¬ 
ture  (  Ts/J  ). 

Since  dissociation  can  proceed  from  any  vibrational  level,  the  assumption 
of  vibrational  equilibrium  imrnediatdy  behind  the  ahock  leads  to  a  greater 
population  of  higher  energy  levels,  and  to  a  higher  rate  of  dissociation.  For 
the  flow  conditions  where  vibrational  nonequilibrium  ie  important,  only  the 
lowest  vihratianal  levels  of  the  molecules  are  populated,  so  that  dissociation 
'.'.ill  proceed  more  slowly  than  if  vibrational  equilibrium  were  established.  In 
se  regimes,  coupling  the  dissociation  rat©  Constant  to  the  degree  of  vibrational 
excitation  is  a  more  accurate  description  of  the  relaxation  phenomenon.  This 
slows  down  the  dissociatioi  process  and  a  higher  temperature  is  maintained 
during  the  period  of  vibrational  relaxation  than  if  the  vibration  and  dissociation 
proceeded  independently. 

Discussion  of  Program  and  Results 

The  normal  shock  program  has  been  completed  for  the  IBM- 704  computer, 
and  has  been  used  to  investigate  relaxation  phenomena  in  a  variety  of'gases.  A 
capacity  of  20  species  and  40  chemical  reaction#  has  been  included  in  the  program. 

In  addition,  the  thermodynamic  description  of  each  species  has  been  improved 
over  that  discussed  in  Reference  1.  Higher  electronic  states  are  now  included, 
with  a  capacity  of  eight  electronic  levels  for  each  species,  adding  considerably 
to  the  accuracy  of  the  solution.  A  variable  interval  is  employed  in  the 

step-by-step  integration  procedure  away  from  the  shock.  In  this  process,  the 
interval  is  doubled  after  a  specified  number  of  successful  steps  (provided 
ceitain  tests  are  passed,  i.  e.  per  cent  cliange  in  T,  change  in  species  concentration, 
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ett.  are  within  given  bounds)  until  a  maxirrium  specified  interval  size  is  reached. 
If  changes  in  an  interval  are  too  large,  the  interval  size  is  decreased  until  all 
tests  are  passed.  Small  intervals  are  then  used  in  regions  of  rapidly  changing 
variables  (i.  e.  near  the  shock  front),  while  progressively  larger  intervals  are 
used  as  the  variables  asymptotically  approach  equilibrium  values. 

The  results  of  a  typical  normal  shock  computation  are  shown  in  Figures 
10  through  12.  For  this  particular  case,  vibrational  equilibrium  immediately 
behind  the  shock  was  assumed.  The  initial  conditions  of  p  =  1  mm  Hg  and  shock 
Mach  number  of  12.  3  are  from  a  typical  shock  tube  test.  Figure  10  shows  the 
species  distribution  as  a  function  of  distance  behind  the  shock  wave.  Seven 
species  (O^.  N^,  O,  N.  NO.  NO+.  e')  were  taken  to  define  the  air  kinetic  model 
and  the  following  six  reactions  were  used:  (  denotes  the  forward  rate 

c  oo 

at  vibrational  equilibrium,  in  cm^/ mole  -  sec  ) : 


1.  O  +  M  •  20  +  M 


2.  +  M  2N  +  M 


3 .  NO  +  M  N  +  O  +  M 


O  +  N. 


N  +  O 


2  ^ 


NO  +  N 


NO  +  O 


6.  N  +  O  NO"^  e‘ 


Ap  =  1.  1  X  10^^  (T) 


25  ,.,,-2.  5 


exp/ 


=  6.  3  xlO^^(T)'^  exp^ilii^ 

"  I0"'(T)-*-  5expDl5iL_ 

w 


Ac  -5.5x10  (T)^’  ^  exp 

■Too 


-1x10  ^'PjO.S 
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The  formation  of  atomic  species  (O  and  N>  i,  seen,  with  the  subsequent  decrease 
in  concentration  of  the  original  molecular  specie.  and  N^)  as  the  distaivce 
from  the  shock  wave  increases.  The  overshoot  in  NO  concentration,  which  is 
typical  of  chemical  relaxation  zones  in  air,  reaches  a  maximum  approximately 
one  centimeter  behind  the  shock.  A  positive  gradient  in  electron  concentration 
exists  until  1.  5  cm  behind  the  shock,  where  the  maximum  electron  concentration 
is  reached.  At  this  point,  there  is  a  20%  overshoot  in  the  electron  concentration 
over  the  final  equilibrium  value.  The  magnitude  of  this  electron  overshoot  is 
dependent  on  shock  Mach  number  and  density,  and  is  as  high  as  a  factor  of  two 
at  the  higher  shock  velocities.  It  should  be  noted  that  equilibrium  concentrations 

are  not  reached  until  some  20  cm  behind  the  shock,  indicating  a  long  relaxation 
zone. 

It  IS  of  interest  to  note  the  role  each  reaction  plays  in  this  relaxation 
process.  An  insight  to  this  can  be  seen  from  Figure  11,  where  the  ratio  of 
total  forward  rates  to  reverse  rates  are  given  as  a  function  of  distance  behind 
the  shock.  Immediately  behind  the  shock,  this  ratio  is  infinite  for  the  dissocia¬ 
tion  reactions,  since  the  reverse  reaction  is  zero,  and  the  whole  process  proceeds 
in  the  forward  direction.  At  final  equilibrium,  this  ratio  would  be  one.  The 
two-body  NO  exchange  reactions  and  the  ionization  reaction  reaches  local 
equilibrium  much  closer  to  the  shock  than  do  the  slower  dissociation  reactions. 

It  IS  the  exchange  reactions  that  produce  the  NO  overshoot,  since  they  tend 
to  put  the  NO  in  equilibrium  with  the  molecular  O^  and  early  in  the  relaxation 
zone.  A  slight  undershoot  in  the  ratio  is  seen  for  these  reactions  indicating 
that  the  rever.so  rate  Miich  deletes  NO)  is  dominant  in  that  regime,  and  the 
decrease  Irom  the  overshoot  maximum  is  aided.  Figure  12  shows  the  effect 
ot  the  chemical  system  on  the  gasdynamic  properties  of  terup* ratufe  and  fle^isity. 
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bchi«,  the  .h.ck.  tv  Ut^e—tete,  .3  S=.„«tne^  ttensletloe, 
tetatloe,  -tlbt.tien  a»d  eleetronit  e.cilatiah,  i.  «ppto«m.tel,  72ao-«.  The 

te„pa„^,te  d.ctee.e.  through  the  r.layation  t.giou  towards  the 

equilibrium  value  of  4400° K  An 

uvj  fv.  An  increase  of  approximately  50%  U  seen  In  the 

density  ratio  during  ,h.  chemical  proce.,.  Thus,  0,e  early  part  of  the  nonequiU- 

hrtum  rone  consists  of  molecuiar  spec...  at  ehtremeiy  high  iransla.ionai  tempera- 

tures.  For  very  high  aUiiude  re-entry,  thi.  early  part  et  the  relaxation  cone 

may  fill  the  entire  shock  layer,  giving  a  s..h.t,et„llv  different  flow  field  than 

that  computed  from  equilibrium  conditions. 

Condition,  of  vibrational  no^tquilibrlnm  have  also  been  investigated 
using  the  norma,  shock  program.  In  «...  analyaie,  bhe  forward  rate  constant 

at  vibrational  equilibrium  for  the  t  reaction  i  %  \  ■ 

V  reaction  {  )  is  modified  by  the 

vibrational  coupling  factor  S/J  through 


\A  ^ 


ij  can  he  1  or  denotmg  chat  vibra, tonal  oquil.brat.on  w.ll  proceed  coupled 
to  or  tudepondeuily  „f  ...e  chemical  reaCtoh.  Since  the  v.brat.onally  relaxing 
.pec.es  arc  considered  be  a  system  of  hanmon.c  oscillators  in  a  Boltamann 

distr.h„tion  about  a  vibrational  temperature.  Tvj  ,  ,hi,  TvJ  will  start 

at  the  original  a.nb.ent  temperature  unn.ediafely  behind  the  shock.  In  general 

'  't  -  0.  Ivj  «  ,  ,  'ib,,  vibrational  relaxation  can  be  seen  in  p-igure  13 

wh.ch  .hows  the  tenperature  d.stribut.on  behind  a  Mach  1,  shock  wave  into  100 
microns  of  air.  For  this  case,  the  degree  of  vibrational  nonequilibrium  was 
oupled  to  the  dissoclafiou  rcacucrns,  and  the  aame  kinetic  air  model  described 
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„3...  Xhe  <e™pe...„e  .3  .e.e.„i„ea  .„„s,a..„„ 

rotat.„„  and  exCed  .leCronic  level,  is  app.«i„a,ely  19,  000-K  at  Al-  =0 
and  decreases  .awards  the  ..juilibriu*  value  of  6000-K.  Both  the  molecular 
otcygen  and  nitrogen  vibrational  teinperatures  are  300- K  a.  lU.  ,0  (the  original 

nt  value)  and  increase  toward,  the  translational  temperature.  Since  the 

ottygen  equilibrates  more  rapidly  than  the  nitrogen  the  TT, 

^  ''^o^  reaches  local 

equilibrium  sooner  fi.  e  ^  ,  -a. 

**  compared  to  ff  =2.  5  cm  for 

V«a,  ».  Since  «,a  translational  temperature  is  decreasing  a,  this^int,  .he 
vtUrational  temperature  overshoot,  th.  local  equilibrium  value,  and  approaches 
e  n^l  equilibrium  result  i,om  above.  This  phenomenon  i.  more  pronounced 
'^0^  *han  for  since  it  occurs  earlier  in  the  relaxatin 

ic  decreasing  verprapdiiy.  When  To.<T,  <  1  and  di,  h. , ward  rate 

constant  is  reduced,  slowing  down  die  rare  of  lorm.tion  oi  .he  atomic  species 
■owever,  during  the  overshoot  regime,  Vj  >  i  .  .„n 

enhanced,  dince  the  atomic  specie.  O  and  «  are  delayed  early  in  the  relaxation 

region,  the  formation  oi  electron,  i.  al.o  delayed  since  it  i  ,h 

^  ^  ^  *6  the  recombination 

Of  and  N  through  ,he  sixth  reaction  dia.  give,  rise  te  ,h.  electrons.  However, 

the  pom.  of  maaimum  electron  concentration  i.  reiatiy.ly  unaffected  by  this 

coupling.  It  is  the  initial  gradient  in  tin  i  a 

gradient  in  the  electron  concentration  which  is  most 

sensitive  to  the  effect  of  \/*  i 

j  .  In  general  then,  the  vibrational  coupling 

pTienomenon  tends  to  <iecrea«ar 

.crea.e  the  dissociation  proc.s.  early  ii,  the  relaxation 

rone,  keeping  the  tr,„.I„i„„ai  temperature  high. 

Conclusions 

1  .hock  solution  gxves  u  ciear  picture  of  the  tJiermocKemical 
P  enomeiia  occurring  behind  .hock  wave.,  rhi.  .oiutlon  cap  he  used  to  gain  » 
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insight  into  the  chejmistry  of  actual  shock  layers  surrounding  re-entry  vehicles, 
since  it  cloaeiy  tesejnbles  the  stagnation  streanaiine  solution.  However,  near 
the  body,  the  actual  flow  velocity  approaches  «cro  along  the  stagnation  streamline. 
Thus,  the  ratio  of  chemical  time  to  particle  reaicJence  time  becomes  small,  and 
the  process  approaches  equilibrium  much  more  rapidly  near  the  body  than 
indicated  by  the  normal  shock  solution. 

The  normal  shock  program,  then,  yields  a  description  of  the  relaxation 
zone,  showing  the  species  distribution  with  the  typical  NO  and  electron  over¬ 
shoots,  the  decrease  in  temperature,  increase  in  density,  and  the  effect  of  the 
coupling  of  vibration  with  dissociation.  This  normal-shock  solution  is  now 
complete  for  any  multicomponent  gas  mixture  of  atomic  and  diatomic  molecules 
and  will  be  used  to  determine  the  importance  of  various  rate  constants,  to  aid  in 
the  interpretation  of  shock  tube  data  and  give  an  insij^h*!  into  the  nonequilibrium 
shock  layer  in  the  stagnation  region  of  hype rvclocity  vehicles  The  computational 
time  is  very  short  for  this  program,  in  that  approximately  45  seconds  are 
required  to  compute  the  relaxation  zone  for  a  single  diatomic  gas  (i.  e.  pure 
oxygen),  and  about  6  minutes  to  obtain  a  solution  for  air  using  the  air  model 
previously  discussed. 

B.  BO  W-SHOCK  SOLUTION 

General  Discussion 

The  tli«‘ rmodynamic  s  and  chemistry  of  this  prografn  are  identical  to  those 
of  the  normal  shock,  and  the  greater  complexity  arises  fnom  a  more  difficult  flow 
situation.  A  shock-wave  size  and  shape  is  specified  ii»  the  invcjtsc  method^'^, 

4nd  the  body  obtained  from  mass  conservation  cc»nsa<ic* Nations  iu  the  •systern  of 
equations.  Consider  tlie  sketch 
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where  X  ,s  the  distance  along  the  shock  and  ^  is  the  coordinate  normal  to  it. 
The  shock  is  divided  into  ho  velues  ot  iX  a.  shown,  and  the  conditions  imme- 
diately  behind  the  shock  at  each  point  arc  determined  from  the  local  normal- 
shock  cohsiderations  discussed  in  the  previons  section.  Once  these  starting 
values  have  been  obtained,  the  entire  solution  advance,  forward  one  interval 
normal  to  the  shock  wave.  The  solution  cannot  proceed  forward  one  ray  at  a 
time,  since  the  derivatives  in  the  X  direction  parallel  to  the  shock  are  needed. 
These  represent  the  flow  swept  by  a  point  on  any  given  ray,  since  the  ray, 
do  not  follow  streamlines,  but  cross  them.  Accuracy  depend,  on  the  accuracy 

with  which  the  X  derivatives  are  known,  the  sire  of  the  mesh  intervals,  the 
method  of  numerical  integration,  etc. 


Discussion  of  Program  and  Results 

Since  the  last  Semi-Annual  Report  ^  was  completed,  the  hi^h- , peed  core 
storage  unit  has  been  installed  in  the  computer,  and  the  bow  shock  solution  has 
been  reprogrammed  to  take  advantage  of  the  larger  capacity.  At  the  present 
time,  this  solution  has  a  capability  for  10  species  and  ZO  chemical  reaction*, 
l^ith  thermodynamic  descriptions  identical  with  those  of  the  normal  .hock 
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A  variable  Interval  haa  alao  been  inccpora.ed  in  the  bo-  .hock  program, 

eliminating  the  need  tor  an  accurate  knowledge  of  the  initial  stop  to  be  used. 

Subject  to  the  various  tests  included  in  the  solution,  the  program  selects  an 
initial  step  size  that  will  satisfy  all  requirements.  Thus,  small  steps  can  be 
taken  both  near  the  shock  wave  and  near  the  body.  For  the  test  cases  discussed 
in  this  report,  a  three-point  fi,  has  been  used  for  the  derivatives.  Upon 
final  refinement  and  completion  of  the  program,  a  more  accurate  five-  or  seven- 
point  fit  will  be  used.  The  method  of  computing  these  derivatives  is  incorporated 
in  a  separate  program  subroutine  that  can  be  easily  changed.  The  programming 
of  the  solution  is  essentially  completed,  and  a  number  of  test  cases  have  been 
run,  both  with  a  diatomic  gas  (i.  e.  O^)  and  with  air. 

1.  Re-Entry  Case 

A  typical  re-entry  condition  was  chosen  as  the  first  air-test  case 

using  the  same  air  model  previously  discussed,  with  vibrational  equilibrium  ’ 

assumed  immediately  behind  the  shock  wave.  A  velocity  of  Z3,  000  ft/sec  at 

200,  000  ft.  altitude  was  chosen.  Using  the  1959  ARnr  cat  u 

using  me  iyt,9  ARDC  atmosphere,  the  ambient 

conditions  are 


T  =  249.  44°K 


^ .  1782  X  10  grn/cc 


■ft  =  2.1281  i  10-4  ^  gm/mole 

Figure  H  shows  the  nose  region  of  the  re-entry  vehicle,  with  a  catenary  shock 
wave  having  an  approwimate  6-.„ch  .radius  a.  the  centerline  specified.  The  hody 
Shape  .as  computed  from  mass-flow  considerations  is  shown  in  the  shock  coordinate, 
.n  Figure  ,1a.  This  plot  actually  shows  how  the  bod,  piiiis  away  from  the  shock. 

A  slight  wavine«s  can  be  seen  in  the  lioHv  <=Vn,-.ca  -tu 

oody  shape,  with  more  pronounced 

inaccuracies  on  the  end  ray..  This  is  due  the  fact  that  a  linear  fit  is  presentiy 
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used  for  the  JJC  derivatives  on  the  end  rays.  These  outermost  points  should  be 
disregarded  for  the  test  cases.  At  the  present  time,  a  program  stop  occurs  when 
the  first  ray  reaches  the  body  (i.  e.  ray  closest  to  centerline),  and  the  solution  is 
not  carried  into  the  body  to  obtain  a  more  accurate  determination  of  the  body 
shape  on  the  outermost  rays.  Thus,  it  is  not  presently  known  whether  the  slight 
waviness  is  due  to  the  3  point  y,  derivative  fit,  or  to  possible  inaccuracies  in  the 
mass  conservation  criterion  as  the  distance  away  from  the  centerline  increases. 

The  entire  JC  derivative  subroutine  is  now  being  refined,  and  this  refinement 
will  lead  to  a  more  accurate  determination  of  body  shapes.  In  the  results  presented 
here,  a  total  of  24  rays  (i.  e.  A^C  divisions)  were  used  in  the  computation.  As 
shown  in  the  axisymmetric  coordinate  system  in  Figure  14b,  the  body  shape 
closely  resembles  a  circle.  The  shock  standoff  distance  is  seen  to  lie  between 
the  computed  equilibrium  and  frozen  values,  which  were  obtained  using  the 
relation  of  li  and  Geiger  ,  The  present  solution  was  carried  to  approximately 
60”  off  the  centerline,  requiring  a  computing  time  of  one  hour  and  fifteen 
minutes . 

Figure  15  shows  the  temperature  distribution  in  the  shock  layer  along  4 
rays  labeled  in  Figure  14b.  The  temperature  immediately  behind  the  shock  as 
determined  from  translation,  rotation,  vibrational  excitation  and  higher  electronic 
levels  is  approximately  17,  500°  K  along  a  ray  near  the  centerline  ^i.  c.  yi  =  .  06 
IS  about  3°  off  the  centerline).  However,  the  temperature  at^=  0  at  approximately 
‘■>0°  off  the  centerline  (i.  e.  X  =•  86)  is  much  less ,  11,000*K.  The  mast  pro¬ 
nounced  changes  in  the  temperature  occur  near  the  stagnation  region  {  %  =  .  (|6) 

since  the  reactions  progress  rapidly  here  due  to  the  high  temperature  and  density 
and  because  the  local  flow  velocity  is  low.  It  is  weSn  that  the  entire  shock  layisr 
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is  in  a  state  of  non.quilibnium,  at  temperatures  greatly  in  excess  of  the  equilibrium 
normal-shock  result  (  ^  6200®Ky, 

Figures  16.  17,  and  18  show  the  species  distributions  in  the  shock  layer 

around  the  nose  cap  of  the  vehicle.  The  Equilibrium  normal  shock  results  are 

given  a,  a  reference  on  each  figure.  These  figure,  should  be  interpreted  as 

giving  the  specie,  concentration  along  a  ray  from  the  shock  (  ^  =0)  to  the  body, 

i.  e.  walking  along  the  ray  normal  to  the  shock  until  the  body  i.  reached.  The 

electron  cohcan.ration  (Figure  ,4)  show,  an  extreme  gradient  from  the  .bock 

to  the  body,  with  an  overshoot  of  some  70%  over  the  equilibrium  value  for 

the  ray  closest  to  the  centerline.  However,  along  an  outer  ray  (  *.  e  ,  gg)  t,,e 

electron  concentration  is  an  order  of  magnitude  less  than  the  value  near  the 

stagnation  region.  The  electron  concentration  gradient  will  have  a  strong  influence 

not  only  on  ga.eou.  radiation  calculation.,  but  al.o  in  ,h.  determination  of  yug., 
cross  -  sections. 

The  chemical  epecie.  distribution  i,  ahown  in  Figure.  17  and  18,  and  it 
is  interesting  ,o  note  extreme  ov.r.boo,  in  NO  concentration  above  the  equilibri, 
mal  shock  solution  (approximately  two  order,  of  magnitude).  This  overshoot 
occurs  on  all  raye,  denoting  that  the  entire  shock  layer  ha.  an  extremely  high 
concentration  of  nitric  oxide  at  the  high  nonequilibrium  temperature.  di.cu,.,d 
previou.ly.  The  oxygen  and  nitrogen  atomic  .pecie.  are  ..e„  to  approach 
equilibrium  value,  within  U,,  .hock  layer.  However,  the  small  oxygen  molecular 
specie,  concentration  i.  still  an  order  of  magnitude  above  the  normal  .hock  value 
When  .he  body  I.  reached.  Thus,  for  this  typical  re-entry  case,  the  front  portion 
•hock  hryer  i,  ,n  a  state  of  nonequilibrium  with  large  amounts  of  NO  and 
et  extremely  high  temperatures.  This  will  give  a  con.iderably  differ.,,,  radial  J 
Picture  than  than  computed  from  equilibrium  cond.lions  (i.  e.  T  »  hOPP-K.  wim 
cimount*  of  NO  artl  0^|. 
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2.  Hypersonic  Tunnel  Case 

The  rapid  expansion  of  the  flo„  from  a  high  temperature  reservoir 
in  hypersonic  test  facilities  may  introduce  noneguilihrium  phenomena  in  the 
expansion  nozzle.  These  effects  will,  in  turn,  have  an  influence  on  the  shock 
layer  around  a  model  in  the  expanded  flow.  Thus,  it  is  necessary  to  obtain 
a  coinplete  flow  solution  for  the  tunnel-model  configuration.  Using  conditions 
similar  to  those  encountered  in  hypersonic  tunnel  investigations  (i.  e.  reservoir 
pressure-  ,0(,atm.  ,  reservoir  tenrpera.ure  =  dOOO'K).  the  finite-rate  expansion 
was  computed  using  a  program  recently  devesoped  at  Cornell^.  The  resultant 
noneouilibrium  flow  condition,  were  used  as  inputs  to  the  bow  shock  program. 

The  results  chosen  test  eection  conditionc  were: 

CoJ  =  .  158  woles/orig.  mole 
^  ^  atnt  =  i.5x  10*^® 

18.i90  tO;j?=.084 


^  -  4.  2  X  10^  cm/sec. 


M  =  .  689 

[ho]  =  .  06 

=  8.6x10 


-3 


Thus,  the  ambient  condition,  CDnsisied  of  large  amounts  of  O  and  NO  indicating 
freezing  of  these  spec,.,  early  the  eapans.on  process.  Again,  a  catenary 
shock  shape  was  assumed,  hav.ng  a  r.diu.  a.  the  centerl.ne  of  approximately 
6  inches.  This  would  presunre.  then,  a  s.x-inch  nos.  zad.us  model  .n  the  nozzle 

test  section.  A,  in  the  pree.ous  test  case,  th,  mod,!  shape  closely  resembled 

a  circle. 

Figure  19  shows  the  tempo  rati, te  di.trihution  in  th.  shock  !aye.-.  There 
i-s  a  marked  difference  in  thi.  distribution  from  that  calculated  for  the  20r^.aor> 
le-entry  case.  riu«  is  due  tv  tlie  low  density  ir  th*  t«st  seciinn  U-  f-.  typicail 
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Of  app„.ta.,e„  250,»00  f..a,.it„de)  logger  relaxation  times,  and  to 

the  presence  o,  large  amounts  of  O  and  NO  in  the  flow.  For  the  innermost  ray 

(>^  =  .  12)  the  temperature  stay,  practically  constant  and  does  not  start  to  decrease 

toward  an  eguilihrium  value  until  the  vicinity  of  the  body.  Along  the  outer  rays. 

the  temperature  increases  slightly,  similar  to  an  isentropic  compression  for 

the  Ideal  gas  case,  before  the  chemistry  begins  to  decrease  it.  For  the 

outermost  ray.  this  latter  effect  is  unnoticeable  up  to  the  point  where  the  body 
is  reached. 

This  pronounced  effect  is  illustrated  in  Figures  20.21.  and  22.  which 
Show  the  species  distribution  in  the  shock  layer.  The  ambient  oxygen  atom 
concentration,  which  is  about  a  factor  of  three  below  the  final  equilibrium  normal 
Shock  value  exists  along  all  rays  well  into  the  shock  layer,  before  an  increase  in 
this  concentration  occurs,  first  along  the  ray  nearest  the  centerline.  However. 

the  body  is  reached  well  before  the  equilibrium  value  is  attained.  The  initial  NO 

concentration  fills  the  entire  shock  laxrfsf  u  • 

ntire  shock  layer,  and  begins  to  increase  in  the  typical 

overshoot  manner.  It  is  seen  that  the  body  is  reached  well  before  the  maximum 

.u  this  overshoot  occurs,  leaving  the  shock  layer  with  an  NO  concentratton  some 

two  orders  of  magnitude  greater  o  ’I't,  • 

g  eater  than  the  equilibrium  normal  shock  value.  The 

IV  atom  concentration  nf*  ^ 

atton  resemble,  the  previous  ease,  since  the  ambient  value  of 

N  is  practically  cero  (i.e.  moles/orig.  mole).  Figure  22  shows  the 

electron  distribution  around  the  mod.l  nose  cap.  It  is  seen  that  the  body  is 

reached  we.l  before  the  overshoot  ft,  electron  concentration  wo.ild  occur.  The 

electron  concen, ration  a.  the  body  is  approx.mately  on.  order  of  magnitude  below 

.be  equilibrium  normal  shock  value.  Thus  the  enUre  shoe*  layer  sarroWin* 

the  nose  cap  is  in  a  •tate  of  «xtireme  nc»ttrquiUbrii«n.  at  r«tu»»,  m,*ch  higV-r 
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than  equilibrium  values.  None  of  the  chemical  species  have  '  :hed  their  equili¬ 
brium  values  at  the  point  where  the  model  surface  is  reached.  The  r  lecular 
oxygen  concentration,  for  example,  is  some  three  orders  of  magnitude  greater 
than  the  equilibrium  normal  shock  solution  in  the  entire  gaseous  envelope.  From 
these  computations,  the  impoftance  of  obtaining  realistic  flow  solutions,  both 
in  the  expansion  nozzle  and  araund  Ji  jnodei.  Is  apparent  for  any  gaseous  radiation 
experiments. 

It  is  interesting  to  note  the  effect  of  noiiequilibriuin  species  in  front  of 
the  shock  on  the  electron  formation  behind  it.  Nornnal  shock  solutions  were 
carried  out  for  the  test  section  conditions,  and  these  results  are  shown  in 
Figure  23.  Curve  A  indicates  the  electron  concentration  behind  a  normal  shock 
in  the  actual  finite  rate  nozzle  test  section  (i.  e.  withO,  N,  NO,  O,,  N  ,  Ar 
in  front  of  shock).  Curve  Cj^  has  the  same  initial  conditipns  as  Curve  A,  with 
the  exception  that  e'  were  included  in  the  finite  rate  expansion.  In  front  of  the 
shock,  a  frozen  concentration  of  b  x  10  electrons/ cc  was  computed.  It  is 
seen  that  the  electrons  start  at  the  ambient  value,  and  asymptotically  approach 
Curve  A  as  the  distance  behind  the  ahock  increases.  The  initial  conditions  for 
Curve  C  consisted  of  the  same  flow  velocity,  temperature  and  pressure  as  the 
preceding  curves,  but  only  the  original  air  composition  was  assumed  (i.  e.  O^, 

N^,,  Ar).  In  this  way,  only  the  effect  of  ambient  nonequilibrium  species  on 
electron  formation  could  be  determined.  The  production  of  electrons  is  delayed 
considerably,  the  concentration  being  some  two  orders  of  magnitude  lower 
at  CO  1  re  spending  distances  behind  the  shock,  for  this  case  as  compared  with 
Curve  A.  Since  electrons  are  formed  through  the  reaction  requiring  the  recom¬ 
bination  of  an  oxygon  atom  with  a  nitrogen  atom,  the  availability  of  these  species 
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immediately  behind  the  shock  enables  the  electrons  to  be  formed  at  a  much  faste. 
rate  (Curve  A). 


C.  SUMMARY  -  ANALYSIS  PHASE 

The  fimte-rate,  normaUshock  program  for  a  multicomponent  gas  ha. 
heen  completed  for  the  IBM-704  computer,  and  has  proven  to  be  an  e.tremely 
helpful  too,  in  ,he  .valuation  of  relaxation  zones  behind  the  shock  waves.  A  great 
variety  of  «...  hav.  been  run  .o  determine  the  effects  of  rate  constant.,  shock 
velocity,  density,  ..c,  on  the  species  formation.  This  program  will  continue 
to  be  used  a.  a  .tasting  point  tor  the  more  complex  bow  shock  solution. 

Th.  general  blunt-body  problem  (i.  e.  the  inviscid,  no„e,jaiiibrium  flow 
over  a  blunt  axi.ymm.tric  body)  ha,  been  programmed,  and  subject  to  fi„l 
refinement.,  .ill  be  completed  in  Ute  near  future.  A  number  of  test  ca...  have 
been  di.cu...d,  .bowing  ,h.  workability  of  th.  program  for  both  re-entry  and 
hyper.onic-tunnel  ...t  configuration.,  A  higher  order  fit  for  ,h,  *  derivative, 
is  being  incorporated  i„m  the  program  to  el.minate  tn.ccoracie.  sn  the  computed 
body  shape.,  b.  addition,  more  eificien,  numerical  tecbnioue.  are  being 
mcorporated  in  th.  program  in  order  to  decrease  the  computing  time. 
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appendix  -  OETERMimXION  OF  THE  EFFECTS  OF 

INTERNAL  REFLECTION  IN  THE  SHOCK  TUBE 


..  .e  .led  .  .„ee..«. 

.K„eL  .Be  I.  Ke.,,.  cBe...,.,.,eTosI.eB  dc.el. 

ph=.pBo,o„,  ^ 

feflection  coefficient  greater  than  .  5  . 

it  has  been  mentioned  in  the  report  that  th.  . 

*  ««Cctive  aperture  into  the 

-OCR  .Be  U  e 

w.^w.  B,c.u..  «..  ..la.  BI  u.. 

ter  of  the  tube,  ,t  behove,  like  ,  cyli„a,ic»l  l.„s.  Figure  9  is  a  sketch 

o  t  .  cros.  eectron  of  the  .hock  tub.  ,,  the  windoo.,  showing  the  accepted  field 

O  view.  Opposite  thi.  window  1,  ..odter  „i,dow  poet,  which,  for  these  studies 
was  plugged  Wid.  .  window  which  ...  ..pp„.,„ 

re.io  doped  Wid.  cebOh  black.  Th„.,  a.e  r.fiec.io,  coefficien.  Of  this  part  of 

•  .  wan  was  decaeod  io  dta.  ei  ih,  .ir-..pphir.  interface.  It  was  measured 

o  c  0.  05.  I,  shouid  he  .entioued  .ha.  for  eea.oo.  of  spectroscopic  cleanliness 
CO  s.uiple  coa„„g  of  pai„t.  p......  .,p., 

»  ecrcase  ,he  wall  reflectanc.  Fro.  tbe  figur.  .,  „„  ,, 

o  .  e  shock  tube  oppos.te  the  wtndow  is  ™w  a  composite  of  plated  „e.al  and 
the  osapphire  window. 

Shock  tube  runs  were  ntade  such  that  the  condi.ioos  of  ih.  gas  within  tb. 

tubf>  were  identical.  fn  one  rH<=«  c  n  • 

"  w.re  us.d  ,1  .  ,g  „„ 

and  m  the  second  the  .lit  was  masked  to  1  x  1  ntb,  Witt, 

iTthn.  With  an  aperture  of  fhi«  size, 
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as  to  limit  the  field 


the  solid  angle  subtended  by  the  spectrometer  was  such 

of  view  on  the  other  side  of  the  tube  to  an  area  completely  bounded  by  the 

sapphire  window.  This  is  also  shown  in  Pi.nre  Tor  .hose  eaperi„on.s 

wh.ch  were  performed  a.  .  9  microns,  a  Dumon.  K-1430  infrared-sensi.ive 

Phoiomnltiplier  was  nsed  to  obtain  a  better  sig„al-to-noise  ratio  than  could 
be  obtained  with  the  PbSe  cells. 

radiation  from  the  gas  was  the  same  in  the  experiments, 
the  ratio  of  the  two  recorded  intensities  were  directly  related  to  the  averaged 
s  within  the  shock  tube.  The  effects  of  varying  angles  of  incidence 
-side  the  tube,  as  well  as  multiple  reflections  are  thus  automatically  compensa¬ 
ted.  This  average  value  of  the  reflection  coefficient  from  the  walls  was  measured 
to  be  0.  4  +  0.  05. 
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